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We continue the study of covariant power-counting renormalizable gravity constrained by scalar 
Lagrange multiplier. Lorentz symmetry breaking is investigated in such a theory in comparison with 
the one in ghost condensation model. Covariant power-counting renormalizable vector gravity which 
is invariant under U(l) gauge symmetry is proposed. Several forms of vector Lagrange multiplier 
in this theory are discussed. It is shown that covariant scalar /vector gravity under consideration 
may have power-law or de Sitter accelerating cosmological solution corresponding to inflationary 
era. Simplest black hole solution is obtained and dispersion relations for graviton are presented. 
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INTRODUCTION 



The building of the consistent and satisfactory renormalizable theory of quantum gravity is fundamental challenge 
for XXI century theoretical physics. Despite the number of attempts (for the introduction, see book [1]) this problem is 
still far from being solved. The idea proposed in ref. [2| for renormalizable quantum gravity is to modify the ultraviolet 
behavior of the graviton propagator in the Lorentz non-invariant way as 1/ |fe| , where k is the spatial momenta 

TT^. and z could be 2, 3 or larger integers. They are denned by the scaling properties of space-time coordinates (x,t) as 

(— i follows, 

x — > bx , t -^ b z t . (1) 

When z — 3, the theory seems to be power-counting UV renormalizable. In the construction of such theory, one 
introduces the terms breaking the Lorentz invariance explicitly (or more precisely, breaking full diffeomorphism 
invariance) by treating the temporal coordinate and the spatial coordinates in a different way. Such model has the 
diffeomorphism invariance with respect only to the time coordinate t and spatial coordinates x transformations: 

^: 8x*=C{t,x), St = f(t). (2) 

Here C(ii, x ) an d /(£) are arbitrary functions. 

In ref. [3j , Hof ava-like power-counting renormalizable gravity with full diffeomorphism invariance has been proposed. 
When we consider the perturbations from the fiat background, which has Lorentz invariance, the Lorentz invariance 
of the propagator is dynamically broken by the non-standard coupling with a perfect fluid. The obtained propagator 
behaves as \/k 2z with z = 2, 3, • • ■ in the ultraviolet region and the model seems to be perturbatively power-counting 
(super-)renormalizablc if z > 3. The price for such renormalizability is the presence of mysterious fluid. Recently 
[J] , this model has been rewritten in covariant form in terms of the scalar field constrained by the specific Lagrange 
multiplier. Hence, the model of covariant and power-counting renormalizable field theory of gravity has been proposed. 
The effective fluid is induced by imposing a constraint for a scalar field by using the Lagrange multiplier field [5j . Due 
to the constraint, the scalar field is not dynamical and even in the high energy region, one can obtain a non-relativistic 
fluid. 

In the present paper, we study the properties of covariant power-counting renormalizable gravities and propose 
new versions of such theories. The emergence of the early-time acceleration which may modify the inflationary era is 
demonstrated. The paper is organized as follows. In the next section we briefly review the covariant power-counting 
renormalizable gravity introduced in ref. [3, |j| • The Lagrange multiplier constraint and flat background expansion 
are discussed. Section three is devoted to the detailed study of Lorentz symmetry breaking in such a theory. The 
comparison with Lorentz symmetry breaking in ghost condensation model is done. In section four the covariant 
power-counting renormalizable vector gravity is proposed. Its different versions are presented and symmetries as 



Also at Tomsk State Pedagogical University 



well as different limits are investigated. Section five is devoted to the construction of accelerating early-time FRW 
cosmology. It is shown that time-dependent accelerating universe or de Sitter universe occurs. The emergence of 
finite-time future singularity and its avoidance via further modification of the theory under discussion is investigated. 
Some outlook is given in the last section. In the Appendix A the scalar-tensor presentation of covariant gravity is 
given where higher derivative gravitational terms are absorbed by additional scalars. Appendix B is devoted to the 
discussion of black hole solution and dispersion relation for graviton. 

II. A MODEL OF COVARIANT AND POWER-COUNTING RENORMALIZABLE FIELD THEORY OF 

GRAVITY 

In this section we review the construction of power-counting renormalizable, covariant gravity introduced in refs.[3j, 
|4J . The action of the model is given by 

S 2n+2 = J d^x^f—g {A - a{{d^d v ^^ v + 2[/ V p V p ) n (d^d^R^ + U R)f 



-A ( ~c^<W + Uo)\ . (3) 



for z = 2n + 2 model (n = 0, 1, 2, • • • ), and 

S 2n+3 = I <£xj—g | A _ a {(a^^Vv + 2U VV p ) n (d^d^R^ + U R)} 

x { (cP WV M V V + 2U V V,)" +1 (d^d^R^ +U Q R)}-\ (^d^d^ + uM . (4) 

for z — 2n + 3 model (n = 0, 1, 2, • • •) (compare with [y]). 

As the generalization, the cosmological constant A term may be added to the action ([3]) and (|4|) as 

S = S 2n +2 - y^2 / d 4x V-9 7 or S = S 2n +3~ 2~2 / ^ x ^f-9- ( 5 ) 

The cosmological constant term does not change the ultraviolet structure of the graviton propagator and therefore 
the power-counting renormalizability. Furthermore, in addition to the cosmological term, the kinetic scalar field term 
may be added 

S = S 2n+2 +Jd i x^g-(-^-±d^cf ) Y or S = S 2n+3 + J d A x^7j (- A _ Id^O^Y (6) 
Note, however, if we shift the Lagrange multiplier field as A — > A — 1, the action ([6]) is transformed as 



S = S 2n+2 + / d*Xyf=9 "TTT + ^O , or S = S 2n+3 + / d 4 xy^g~ [-— + U . (7) 



2~^ + Uo )' ° r S = S ^+3 + j d^^^i^- — 
Then as a result, the cosmological constant is shifted by 

A -> A - 2k 2 U . (8) 

Hence, the ultraviolet properties do not change if the kinetic scalar field term is included. 

Note that the above actions © and (QJ (and ([5]) and ([5])) are totally diffeomorphism invariant as they are given in 
terms of the local fields. In the actions ((3]) and (J4J ( ([5]) and (O), A is the Lagrange multiplier field, which gives a 
constraint 



-d^d^ + u = o, (9) 



that is, the vector (d^cj)) is time-like. At least locally, one can choose the direction of time to be parallel to {d^4>). 
Then Eq. © has the following form: 



Then the spacial region becomes a hypersurface where is a constant since the hypersurface is orthogonal to the 
vector (da<f)- 

On the other hand, by the variation of <f>, for example for z — 2 (n — 0) case in (j3|), we find 

= 4a0" {d v 4>R^ v {d p 4>d a ^R pa + U Q R)} + 0" {Xd^) . (11) 

For z > 3 (n > 1 case in Q or n > case in |@J), one gets rather complicated equation. 

Let us consider the perturbation from the flat background g^ v = r\ pv + h^ v . Then the curvatures have the following 
form: 

iV = \ [dn&Kp + d»d p K P - d p d p h^ - d p d v (v pa h pa )} , r = d»d v h pv - d p d p (ifhpa) ■ (12) 



The following gauge condition is chosen: 



Then the curvatures (fT"3l) look as 



ht = hi = ha = . (13) 



Rtt = -\d 2 t (6 ij hij) , Rij = \ {did% k + dj&hu, + d 2 t h l2 - d k d k h ZJ } , 

R = tf&h^ + d 2 t (tfihij) - d k d k {SPhij) , (14) 



and wc find 



d^d^R^ + UoR = U {d i d j h ij -d k d k (6 i 'h ij )} , 
d' i <t>d u (l>V li V v + 2U V p V p - 2U d k d k . (15) 

In the ultraviolet region, where k is large, the propagator behaves as 1/ \k\ for z = 2 (n = 0) case in ([3]) and therefore 
the ultraviolet behavior is rendered. For z = 3 (n — 0) case in Q, the propagator behaves as 1/ |fc| and therefore 
the model becomes power-counting renormalizable. For z — 2n + 2 (n > 1) case in ([3]) or z = 2n + 3 (n > 1) case in 
0, the model becomes power-counting super-renormalizable. The dispersion relation of the graviton is then given by 

u> = acok z , (16) 

in the high energy region. Here Co is a constant, u is the angular frequency corresponding to the energy and k is the 
wave number corresponding to momentum. If a < 0, the dispersion relation becomes inconsistent and therefore a 
should be positive. 

In this section, we only considered the perturbation from the flat background which has the Lorentz symmetry. 
In principle, any manifold is fiat in a microscopic limit (which is a part of the definition of "manifold") and the 
renormalizablity is the problem in such a short distance region. In the microscopic limit, all the manifold can be 
regarded as fiat. Thus, the renormalizability should be discussed only in the flat background. 

Note that any gravity theory does not respect global Lorentz symmetry since the theory is formulated in curved 
space-time. Only when we consider the perturbation from the flat background, which has a global Lorentz symmetry, 
we can discuss about the Lorentz symmetry. Then what we need is a model which admits a solution where the space- 
time itself is flat and therefore Lorentz invariant but a vector and/or tensor has a non-trivial value, which breaks the 
Lorentz symmetry. In this section, (d p 4>) is used as a vector. More general case will be discussed below. 

The above construction may be generalized for D dimensions. If z = D — 1, the model becomes power-counting 
renormalizable and if z > D — 1, the model becomes super-renormalizable. For example, for D = 11, which may 
correspond to the M-theory, if n — 4 in Q is chosen, the power-counting renormalizable theory emerges. 

III. STRUCTURE OF THE LORENTZ SYMMETRY BREAKING 

Let us study Lorentz symmetry breaking for the actions ([3]) and (|4]). The actions are manifestly invariant under 
the full diffeomorphism invariance, which is a local symmetry. The action has a shift symmetry 

^^0 + 0o- (17) 

Here &o is a constant. 



One can confirm that the actions admit a flat vacuum solution. Indeed, field equations are: 

o = -^ U„u - \g^ + c^ ghcr - \d^dA + \ 9ixv Qd p w + Uo^j . (is) 

Here G°^ comes from the higher derivative term (the second term) in the actions ^ and (£[]). When we assume 
the flat vacuum solution, Eq.© given by the variation over A has a form (|10[) . Since for the flat vacuum solution all 
the curvatures vanish, Eq.(fT8|) reduces to 

= \d fl <j>d v <t>, (19) 

whose solution is A = since d^<p 7^ due to the constraint equation @ (if we choose the coordinate system properly, 
we have dt4> — V%Uo arj d di4> — 0). Hence, the actions ([3]), and ((4]) admit the flat vacuum solution with A = 0. 

We now consider the perturbation from the flat vacuum solution g^ = -q^ + h^ as in (|12j) . Then the actions (|3|) 
and (U]) obviously enjoy the Poincare symmetry, that is, the actions are invariant under the Lorentz transformation 
and the translations. In the flat vacuum background, a solution of (| 10[) is given by 

<f) = y/2U^t . (20) 

Since the solution depends on the time coordinate i, the solution spontaneously breaks the symmetry under the time 
translation: 

t^t + to (21) 

Here to is a constant. The solution (|2"0")) also breaks the shift symmetry in (fTT|) . We should note that the diagonal 
symmetry of the time translation (|2"Tj) and the shift symmetry (|T7|) is not broken. In fact, if we choose 4>q in (fT7|) as 

cf>0 = -y/2Uo~t , (22) 

under the simultaneous transformation, the solution (|20j) is invariant. The diagonal symmetry effectively plays the 
role of the time translation and the fiat vacuum solution is effectively invariant under the time translation (for rigorous 
discussion of Lorentz symmetry breaking in effective theories, see Q)- This structure is almost identical with that in 
the ghost condensation model Q, which is a variation of the k- inflation Q or k-essence [lCJ models. The ghost sector 
of the ghost condensation model is given by 

S gh o S t = J d 4 xV=g~P (g^d^d^) . (23) 

Here P (g^d^dvrfi) is an appropriate function of g^d^d^ip. By the variation of (j>, one obtains 

= d^ {^f~gg^d^P' (g^d^d^)) . (24) 

In the flat background g^ w — rj^, a solution of (J24I) is given by 

4> = fot. (25) 

Here /o is a constant. The solution (|2"5|) is identical with that of ([2"U| if we identify /o = \/2Uo- Since the vector d^cj) 
has a non-vanishing value, the Lorentz symmetry is broken in both of the ghost condensation model and our models 
© and (HJ). There is, however, a big difference between the ghost condensation model and our models. When we 
consider the fluctuation from the solution (j2"5j) in the ghost condensation model <fr — f t + <50, the fluctuation has a 
propagating mode. In our models, however, if we denote the fluctuation from the solution (|2"0"|) as <f> = ^/2Uot + 5<f), 
Eq. © shows that 



whose solution is given by 



at =0 ' (26) 



H = f («0 ■ (27) 



Here / (x) is an arbitrary function of the spatial coordinate x. Since S(j> does not depend on the time coordinate, 
there is no oscillating mode. Due to the existence of the oscillating mode in the ghost condensation model, the 



Lorentz symmetry restores in the ultraviolet region and the breakdown occurs only in the infrared region. On the 
other hand, since there is no propagating mode of the scalar field <f> in our model, the Lorentz symmetry breaking 
occurs even in the ultraviolet region. Since we like to modify the propagator of the graviton in the ultraviolet region, 
the symmetry breaking should survive even in the ultraviolet region. This situation is very different from that in the 
usual spontaneous symmetry breaking. 

In the ghost condensation model or our models, the Lorentz symmetry breaking is spontaneous. The usual U(l) 
Higgs model, whose potential is given by 

^Higgs^-^^+y^) 2 , (28) 

has a global U(l) symmetry, which is the invariance under the transformation 

(j)->e ie °<j), (29) 



with a constant real parameter 9q. In (|28|) . is a complex scalar field and m and A are positive parameters. The 
minimum of the potential is given by 

<t> = -^- (30) 

Here ip is a constant phase. The value of tp can be arbitrary. If one chooses specific value of tp, the value of p is 
changed under the t/(l) transformation (|29[) as ip —></? + 9, and therefore the ground state is not invariant under 
the U(l) transformation and the U(l) symmetry breaks spontaneously. One can always choose the real axis of the 
complex 0-plane to be parallel with the value of (j> in the ground state so that tp = 0. 

In our model, the constraint equation © shows that the value of the vector (d^cj)) is located on the hyperboloid 
defined by 

- Z% = t 2 - x 2 = 2U . (31) 

The value of the vector (d^cf)) changes on the hyperboloid under the Lorentz transformation. If we choose a value of 
(dfj,(f)), the Lorentz symmetry is broken spontaneously. After that we can always choose the time axis to be parallel 
to the vector (<9 p </>)- 

Thus, there have been proposed several models, for example, ghost condensation one, where the Lorentz symmetry 
is broken spontaneously but only in the infrared region. We need the Lorentz symmetry breaking in the ultraviolet 
region, where the renormalizability issue becomes a problem. At least, the authors do not know such a model except 
the one in this paper, where the Lorentz symmetry is spontaneously broken even in the ultraviolet region. 

Note that the Lorentz symmetry breaking is not directly related with the renormalizablity. Only if the symmetry 
breaking improves the ultraviolet behavior of the propagator as in this paper, the renormalizability properties are 
changed. For vector theory with spontaneously broken Lorentz symmetry, as in the next section, the explicit coupling 
improves the renormalization properties. 

IV. COVARIANT VECTOR GRAVITY 

The spontaneous Lorentz symmetry breaking occurs when the quantity like vector or tensor, which is generally not 
invariant under the Lorentz transformation, has non-trivial vacuum expectation value. An exception occurs in case 
of the rank 2 symmetric tensor B^ v : if the vacuum expectation value of -B M „ is proportional to r\^ v '. B^ = cr\^ v with 
a constant c, the Lorentz symmetry does not break since rj^ v is invariant under the Lorentz transformation. Then 
instead of the actions © and ((H) , by using a vector field A^ , we may consider 

S 2n +2,A = f d^x^—g | A _ a {(AWVuVv + 2U VV p ) n (AM% + U R)} 2 



-A ( -A^A" + U ) \, (32) 



for z = 2n + 2 (n = 0, 1, 2, • • • ) and 

S 2n+S , A = / d^xj—g | A _ a {(A"A W V^V„ + 2U V p V p ) n (A»A»R^ + U R)} 

x {(AM"V„V„ + 2f/oVV p )" +1 {A^A"R yv + U R)} - A QiM" + U ) } , (33) 



for z = In + 3 (n = 0, 1, 2, • • • ). A possible problem is that A M can take a random value as long as it is time- like at 
different points of the space-time. This could be compared with the vector (d M </>), which is almost parallel even at 
different points of the space-time. At least if we start with the background where the value of the vector field is parallel 
at different points in the space-time, which is also a solution, the gravity seems to be power-counting renormalizable. 
In fact, the constraint equation given by the variation of the Lagrange multiplier field A is given by 



1 
2 

A solution of (|34|) is given by 



= -ApA" + U . (34) 



A = V2£/o, A=0, (i= 1,2,3), (35) 

which breaks the Lorentz symmetry spontaneously. Then 

AWR^ + UoR = U Q {d i d j h ij -d k d k {5 ij h ij )} , 
A»A v V tl V v + 2U W>V p = 2U d k d k , (36) 

which gives expressions identical with (|15[) and therefore we obtain (super-)renormalizable theories for z = 3 and 
z = 2n + 2. 

In general, one can add the kinetic term for the vector field and require the U(l) gauge symmetry. Then the 
following model may be proposed: 

d^V^Q (A - a {(( A>i - 9A V) ( A " ~ V<P) V M V, + 21/oVVp)" ((A" - 0*V) {A" - d v ip) R, lu + U R)} 2 

~ {d^A u - d v A^) (0"A" - PA") - A Q (A M - d,#>) (A» - d^) + U ) } , (37) 

for z = 2n + 2 (n = 0, 1, 2, • • • ) and 

S = / d 4 *^ [ ^ - a {((A» - d»ip) {A" - d v ip) V M V, + 2[/ V%)" ((A" - #V) (A" - 6T<p) R» v + U R)} 

{((A" - 0*V) (A" - 9» V M V„ + 2f/ V' , V p )" +1 ((A" - d»<p) (A v - d v <p) R^ + U Q R)} 



x 

-^ (3 M A„ - dvAJ {d»A v - d v A») - A Q (A„ - 9 M ^) (A" - &><p) + Uq) } , (38) 

for z = 2n + 3 (n = 0, 1, 2, • • • ). Here e is the gauge coupling constant and the scalar field ip is the Stuckelberg field 
The actions (f3"T|) and (j3"5)) are invariant under the t/(l) gauge transformation 

Ay, -»■ A^ + 9 M e , 93 -> 93 + e . (39) 

Here e is the (local) parameter of the the 17(1) gauge transformation. Especially if one chooses unitary gauge 

cp = . (40) 

the actions (l37l) and (l38t reduce to 



5 = y ( i 4 a ;y^|^-a{(A^^V M V, + 2C/ V' , V p ) n (^A^^ + C/ i?)r 



^2 ( 5 a^ - ^A M ) {d»A v - 9" A") - A (~ 



— (0„A„ - flU„) (0M" - a"A") - A -A^A* + U Q ) } , (41) 



for z = 2n + 2 (n = 0, 1, 2, • • • ) and 

S = Jd i x^l^-a{(A»A»V^„ + 2U VV p ) n (A^A»R lw + UoR)} 



{(A"A"V M V„ + 2[/ VV p ) n+1 (A"A"i^ + U R)} 
^ (d„A u - d„AJ {d»A v - d v A») - A Q. 



-— (fl M A„ - 5„A M ) (d»A» - d v A») - A -ApA* + U ) \ , (42) 



for z = In + 3 (n = 0, 1, 2, • • • ). For the solution (|35|) of the constraint equation ([34]) . the held strength d^A^ — d^A^ 
vanishes and therefore the corresponding energy-momentum tensor vanishes. This shows that the flat space-time is a 
solution of the theory (|37| and ([38]) . We should note, however, there appears a constraint given by the variation of <p, 
which corresponds to the Gauss law constraint in QED. Due to the constraint, the models (|37p and (j3"5)) are different 
from the models (j32|) and ([33]) . 

For the gauge fixed action (J4TJ) and (|42"T) in the ultraviolet region, since we have the solution (|35|) . the propagator 
behaves as 1/ |fe| for z = 2 (n = 0) case in (|4Tj) and therefore the ultraviolet behavior is rendered. For z — 3 (n = 0) 
case in (|42|) . the propagator behaves as 1/ |fc| and therefore the model becomes power-counting renormalizable. For 
z = 2ri + 2 (n > 1) case in (|41|) or z = In + 3 (n > 1) case in (j42|) . the model becomes power-counting super- 
renormalizable. 

One may consider the "weak" coupling limit where e — > 0. Then in order for the action to be finite, the field 
strength should vanish: d^A u — d^A^, which indicates that the gauge field is the pure gauge. Therefore, the gauge 
field can be rewritten as A^ = d^4>. By substituting this expression into (J4TJ) and (j42|) . the actions ([3]) and (j4} are 
re-obtained. On the other hand, one may consider the "strong" coupling limit where e — > oo. In the limit, the kinetic 
term of the gauge field vanishes and the actions in (|3"2")l and (1331) are reprouced. Note that the renormalizability does 
not depend on the magnitude of the coupling. 

Now, the role of constraint (f34|) may be investigated. Instead of using the Lagrange multiplier field, the constraint 



can be realized by considering the following action: 




Here Ao is a constant. In the limit of Ao — > oo, the constraint (f34j) follows. The mass tha of the vector field A^ is 
given by 

m 2 A = Aot/o ■ (44) 

Instead of the limit of Ao — > oo, it could be enough to choose Ao so that rriA could be a cutoff scale for the renor- 
malization. In case of the string theory, the natural cutoff scale is the Planck scale. We should also note that there 
appears an infinite tower of the particle modes whose masses are of the Planck scale order. Then in the string theory, 
the cutoff scale mass appears always. In the original proposal by Hofava [2J, the gravity model was expected to be 
an effective theory from the string theory. Then the cutoff scale mass presence looks quite natural. 

In principle, the constraint (|34[) should be imposed only in the high energy region for the graviton. Then we may 
consider the following term: 

Sx G =-^ I d 4 x^—gR n (\A^ + U ) . (45) 



2 

Here Ag and n are positive constant. When graviton has high energy, R becomes large, and therefore the constraint 
P4p appears in the high energy region. It is interesting to remark that Lagrange multiplier constraint which breaks 
Lorentz invariance in the analogy with massive gravity (for recent review, see [ll|) with Lorentz-symmetry breaking 
masses maybe also proposed. 

One may also consider the following model: 

S = S 2n+2 + J d A x^F(G) , (46) 

for z — In + 2 (n = 0, 1, 2, • • • ) model and 

S = S2 n+s +[d i xV=gF(G), (47) 

for z = 2n + 3 (n = 0, 1, 2, • • • ) model. Here G is the Gauss-Bonnet invariant defined by 

G = R 2 - AR^R^ + R^ pr7 R^ pa . (48) 

The F(G)-term in (|4"6")l and (|4"T]) does not change the ultraviolet structure of the propagator of the graviton and 
therefore the models (|46[> and (|47[) remain to be power-counting renormalizable. Thus, the construction of number of 
power-counting renormalizable vector gravities is explicitly presented. 



V. ACCELERATING FRW COSMOLOGY 



The gravitational terms different from general relativity in ([3]) and Q are relevant in the high energy region. 
Such terms might affect the inflationary era. In this section, we briefly study FRW cosmology in the theory under 
discussion. In order to obtain the FRW equations, the following form of the metric is assumed: 

v 2 



ds 2 =~e 2b ^dt 2 +a(t) 2 J2 ( da 



i=l,2,3 



and that the scalar field <j> only depends on time. Then the Eq.© has the following form: 



\{t)=^"° 



Hence, one gets 



d^(j)d v (j)R^ + U R = 6U e- 2b H 2 , ^^"^V„ + 2£/ V p V p = -6U e~ 2b Hd t . 



and the actions © and @ have the following form: 

S 2n+2 = Jd 4 xa 3 ^(6H+12H 2 -6bH)-(6Uo) 2n+2 e b {(e- 2b Hd t ) n (H 2 e- 2b )y 



S> 



2n+3 



-A — 



4„„3 



d xa 



2 \dt 
2k 2 



e b U 



QH + 12H 2 - 6bH 



-2 2n+3 • 3 2 a U 2n+2 e b {(e~ 2b Hd t ) n (tfV 26 )} Ue' 2h Hd t ) n+l (H 2 e - 2b )} 



2 V dt 



e b U 



The first FRW equation looks as 

3 



= ^H 2 - (6t/ ) + aa~ 3 a 3 (D n (H 2 )Y - A{-l) n H 2 D n (a 3 D n (H 2 )) 



n-l 



-4^ (-I) fc (D n - k (H 2 )) (D k (a 3 D n (H 2 ))) 



k=0 



— 2XUq — /^matter j 



for (1521 and as 



= ^H 2 - (6U Q ) 2n+3 aa~ 3 [a 3 (D n (H 2 )) (D n+1 (H 2 )) - 2(-l) n H 2 D n (a 3 D n+1 (H 2 )) 

n-l 

-2(-l) n+1 H 2 D n+1 (a 3 D n (H 2 )) - 2^ (-I) fe (D n ~ k (H 2 )) (D k (a 3 D n+1 (H 2 ))) 



k=0 



-2J2(-V k (D n+1 - k (H 2 )) (D k (a 3 D n (H 2 ))) 



fc=0 



2A(7o — Pmatter > 



(49) 

(50) 
(51) 



(52) 



(53) 



(54) 



(55) 



for ([5"3")l . Here /O ma tter is matter energy-density. We also have put 6 = after the variation over b, where the metric 
reduces to the standard FRW metric and the operations of D and D for a scalar ip are defined by 



n - u di f h - d(H<p) 



(56) 



On the other hand, by the variation over a, we get 
f 







2H + 3H' 



2 2n+2 3 2n+1 aU 2n+2 



3 (D n (H 2 )Y + 4 (-!)" a~ 3 ^ (HD n (a 3 D n (H 2 ))) 



dt 



+2E< 



k=l 



d 

dtWdb 



D n-k (jj2^ \ (£,fc-l ^ D n (#2))) 



' Pmattcr j 



(57) 



for (El) and 



= ^(2H + 3H 2 ^+2 2n+3 3 2n+2 aU^ +2 {-3(D n (H 2 ))(D n+1 (H 2 )) 

+2 (-1)" a- 3 -^- (HD n (a 3 D n+1 (H 2 ))) + 2 (-1)" +1 a~ 3 ^- (HD n+1 (a 3 D n (H 2 ))) 

fc=l ^ ' 

+ f>~ 3 | ((I ( Dn - k+1 (# 2 ))) i^ 1 ^° n (# 2 ))))} +h*«. ( 58 ) 

for (j53")l . The simplest case is n = in (|52")l when FRW equations are 

-iff 2 = -108aC/ ( 2 H 4 + 2At/o + p matt e r , (59) 

—^(2H + 3ff 2 ) = 36aC/ 2 (3ff 4 + 4ff 2 ii) + p matter . (60) 

At the early universe where the curvature was large, the contribution from the Einstein term, which corresponds to 
the right-hand side in (|59")l and (|5(I)) . and the contributions from the matter p m attcr and p ma ttcr, could be neglected. 
Then a solution of (p0|) is given by 

*-5- (61) 

which expresses the (power law) accelerated universe expansion corresponding to the one with perfect fluid with 
w = -1/2. Eq.([39]) gives 

A=^. (62) 

One may consider the actions (j3|) and (j4)), which contain a cosmological constant A. Then the equations ([59)) and 
(1601 look as 



"k 2 1 



— ff 2 = -108aC/ 2 J ff 4 + 2At/ + TT + Pmatter, (63) 

2H + 3H 2 ) = 36aU 2 (sH 4 + AH 2 h) - — + p matter . (64) 



In case that the contribution from the matter is neglected, there emerges the de Sitter solution with constant H = Hq, 
Here Hq is given by solving the following algebraic equation: 



A^_A 

K 2 ° 2k 2 



= 108aU 2 H 4 + —Hi - — ^ , (65) 



whose solution is given by 



H ° ~ '^J 2 ^ + V [aU 2 ^) + 216a[/ 2 r 2 " (66) 



Hence, de Sitter universe which may correspond to inflationary era occurs as the explicit solution. The study of 
cosmological perturbations here maybe done in the same way as in Hofava-Lifshitz gravity (see, for instance, ref.[12j|). 
In addition to the cosmological term, one may add the scalar field kinetic term as in ([6]). Then the equations (J59|) 
and (|60|) are modified as 

^H 2 = -108aU 2 H 4 + 2\U Q + U^\ + A + Pmatter , (67) 

~(2H + 3ff 2 ) = 36aU 2 (sH 4 + 4ff 2 ff) + I f^\ - A + Pmatter . (68) 
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If the contribution from the matter is neglected, there appears again de Sitter universe solutions 



if 



Note that A is given by 



72 // 72 * 


\ 2 A-2k 2 [/ 2 
/ + 216aU$K 2 


A > 4k 2 C/ 2 . 





(70) 



x dty _ i 
A -"^T~~2- (71) 

We should note that the shift of the cosmological constant © , corresponding to the shift of the Lagrange multiplier 
field A — >• A — 1 , is consistent if we compare (l69t with (l66l) . 

In the presence of matter fluid with the equation of state (EoS) parameter w, matter energy-density p m atter behaves 
as 

Pmatter ~ O^ 1 ^ . (72) 

Then if w < — 1, the matter energy-density p ma ttcr increases with the expansion of the universe. In the Einstein 
gravity, when w < — 1, the Hubble rate behaves as 



"Einstein ~ ~ ~ ) (?3) 

t s t 

since -ff 2 oc pmatter • Here t s is a constant and therefore there appears a singularity at t — t s , which is called "Big 
Rip" singularity. In case of (|63| . near the singularity, the Hubble rate H becomes large and H A term dominates and 
therefore we have H oc Pmatter, whose solution is given by 



H ~ 3(1 _+" ) . (74) 

Then there still appears a singularity at t = t s . Therefore the higher derivative term does not prevent the singularity 
unlike to the case of modified F(R) gravity where i? 2 -term cures the singularity as it was observed in refe.[13J. 
In order to avoid the singularity, one may add the following term to the actions ([3]) and @: 

S 2n +2 ->■ S 2n +2 + S A , S 2n+3 ->■ S 2n + 3 + S A , S A = -A J <£xj—g (d^d^R^ + U R) m . (75) 

In high energy region, the mass dimension of d xdt becomes —3 — z. Then the operators with the mass dimension 
less than or equal to 3 + z are power-counting renormalizable. From the actions J3J and Q, it is natural to assume 
the mass dimension of the scalar field is — 1. Then the dimension of the operator (d^ l (f>d v (f>R flv + UqR)™ 1 is 2m. Then 
we assume 2m < 3 + z = 5 + 2n for the action ^ or 2m < 6 + 2n for the action Q. As clear from (fTS")) . when we 
expand the term in the powers of hij, the series start with the TO-th power of h%j. Therefore if m > 3, the term does 
not change the propagator and therefore it does not affect the UV structure. In the background (|4T))) , using ([5TJ) one 
gets 

S A = -A f d 4 x (6U ) m e { - 2m+l ^H 2m . (76) 



Then in n = case, Eq. (|59|) becomes 

3 



H 



2 H 2 = -108aU^H 4 + (2m - 1) A {W ) m H 2m + 2XU + p ma ttcr , (77) 



If there occurs a Big Rip type singularity, the second term of the r.h.s. and the matter density p ma ttcr would dominate 
and we would obtain 

- (2m - 1) A (6C/ ) m H 2m + p matter . (78) 
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However, the quantities in the r.h.s. are always positive if Aq is positive, which leads to the inconsistency. Therefore, 
the Big Rip singularity does not occur. When curvature is large, if the first and second terms in the r.h.s. in ([77)1 
dominate, de Sitter space universe occurs. In the de Sitter space, the Hubble rate is given by 

H=( ™*% Y^ (79) 

The above FRW cosmology is also realized for the theory with the actions (|5^|) and ([55]) where the vector field is 
included, The solution of the constraint equation (|M|) is given by, instead of (j3"5)) , 



A Q =S t) ^2U^, A t = 0, (i = 1,2,3). (80) 

By using the solution (JSUJ), we obtain (j54"l) . (1551) . (I571) . and (1551) . again. Therefore, FRW cosmological solutions do not 
change from the cases of covariant gravity with scalars. One may also add the cosmological term as, 

S = S 2n +2,A - TT-, / ^XyJ^g , or S = S 2n +3.A - -pr-9 / d 4 xy/^g. (81) 

2k J 2k J 

Such extra term does not change the ultraviolet structure of the graviton propagator and therefore, power-counting 
renormalizability. The cosmological solutions are again identical to the ones obtained in this section. 

Thus, we demonstrated that covariant power-counting renormalizable gravity naturally leads to accelerating early- 
time expansion which may correspond to inflationary era. However, the accelerating cosmology is the same for scalar 
or vector covariant gravity under discussion. 

VI. DISCUSSION 

In summary, we constructed and investigated the covariant power-counting renormalizable scalar and/or vector 
gravity which is constrained by scalar/vector Lagrange multiplier term. The ultraviolet behavior of the theory is 
improved because of the Lorentz symmetry breaking which occurs due to non-trivial coupling with the effective 
(scalar/vector) fluid. The comparison of Lorentz symmetry breaking in the theory under consideration with the same 
in ghost condensation model is done. The U(l) gauge symmetry structure of the vector covariant gravity is studied 
as well as weak and strong electromagnetic coupling constant limits. As proposed theory pretends to improve the 
ultraviolet behavior of gravity at high energies, the early-time FRW cosmology is discussed. We demonstrate that 
early-time inflation may occur in the theory under consideration. It is described by power-law accelerated Hubble 
rate or by de Sitter universe behavior. It is interesting that power-law accelerating FRW evolution may end up 
at finite-time future singularity. Nevertheless, the additional modification of the covariant action by extra higher 
derivative term which does not destroy the good ultraviolet behavior of the graviton propagator may cure the future 
singularity. Some remarks about dispersion relations compared with general relativity ones are also made. 

The theory under discussion may be considered as some step towards to the construction of renormalizable quantum 
gravity. Indeed, the ultraviolet behavior of such theory is improved in the same sense as the one in Hofava-Lifshitz 
gravity where Lorentz symmetry is broken from the very beginning. In the covariant theory under consideration, 
Lorentz symmetry is broken dynamically as the result of the non-trivial coupling with the effective fluid. In order 
to understand better its renormalization properties, the one-loop renormalization should be done. This is quite 
non-trivial task due to presence of Lagrange multiplier. It will be considered elsewhere. 
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Appendix A: Scalar-tensor representation of covariant action 

We now consider the structure of the higher derivative terms in the actions ([3]) and ((4]) . First we should note that 

d^d^R^ + U R = (d» <&"<!> - ^d^d'A R^ = (r^ - ~g^R\ d^d v <t>, 
d^d^V^Vu + 2C/ V p V p = (d»<t>d"<f> - g^d p ^d p ^) V M V„ = d^d^ (V M V„ - g^V p V p ) . (Al) 

Then in case of n — in (J3]) , by introducing a scalar field rj, one can rewrite the action as 



S->= I d 4 x^gl^ + aL 2 -2T 1 (R^~^g lw R\d^d v cf,\ A ( -<>„< ><)",. + r„ ) \ ■ (A2i 



Note that the term like (i? M „ — ^g^R) d^<\>d v <\> appears as 0(a) correction in superstring theory. For general n in 
([3]), the action may be rewritten in scalar-tensor form with 2n + 1 scalar fields t], £j, and Q (i = 1, 2, • • • , n): 

S 2 = / d 4 :r^ | A + a | 7? 2 _ 2rid^d v (f > (V M V„ - <? M „ VV p ) Ci 

-A Q^<OT + U ) | . (A3) 

Similarly for ([3]), scalar-tensor action with In + 2 scalar fields £, and £, (z = 1, 2, • • • , n + 1) looks as 

S 2 = y"dW^j ^ + a|ciC2 + ^^(Ci-^^(V M V I/ -^VV p )Cm) 

+Cn+i fcn+i - 6v - \g^R\ d^d'A | - a Qd^> + uM . (A4) 

Appendix B: Black hole solution 

The arguments around Eqs. (JT5J) and (|19p demonstrate that the Schwarzschild black hole 

ds 2 = - (l - ¥\ dt 2 + (l - y^j dr 2 + r 2 d0 2 , (Bl) 



and Kerr black hole 



ds\ = Adt 2 -^rdr 2 - Y?dQ 2 - ^-J- (dip - fldt) 

_ £ 2 (A - a 2 sin 2 9) „ 2aMrA , , , 9 9 9 „ . , 

fi = — — — , A = r 2 - 2Mr + a 2 , £ 2 = r 2 + a 2 cos 2 (9 , (B2) 



" E 4 A - 4a 2 M 2 r 2 sin 2 (9 ' " £ 2 



are solutions since the curvatures vanish. In (|B1[) . d£7 2 is the metric of two dimensional sphere with unit radius. 

For the Hof ava gravity and the theories under consideration, the dispersion relation of the graviton is given by (|16[) 
in the high energy region. Then the phase speed v p and the group speed v g are given by 



u 



du> 



v„ = — = acok z 1 , v R = — — = acozk z 1 , (B3) 

k dk 

which become larger and larger when k becomes larger and goes beyond the light speed. This tells that even in (IB1[) . 
the high energy graviton can escape from the horizon. Note that the horizon is null surface and therefore in the usual 
Einstein gravity, particle cannot escape from the horizon since the speed of the particle is always less than or equal 
to the light speed. In our model, however, the speed of the graviton can exceed the light speed and escape from the 
horizon. 
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The dispersion relations in (IB4|) and /c-dependent speed of graviton could give some effects at the early universe. 
At the early universe, the quantum fluctuations with the wave length A = l/2nk become classical ones when the 
wavelength is larger than the horizon radius rn: 

*-55 >,B -!- (B4) 

Here H is the Hubble rate. In case of (|B4j) . since the speed depends on k, the horizon radius depends on k as 

^ = 1 = ^. (B5) 

Then Eq. (|B4|) shows 

2TTc zk z < H . (B6) 

which tells that the fluctuations in the high frequency modes could be changed from the usual Einstein gravity 
2irk < H (we choose the parameter "light speed" to be unity c = 1). The fluctuations could be observed via CMB. 
Thus, in future the possible difference in the high frequency region might be observed what may provide the additional 
test in favor/against of theory under discussion. 
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